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Abstract 

Spin-polarized neutron matter is studied using chiral two- and three-body forces. We focus, in 
particular, on predictions of the energy per particle in ferromagnetic neutron matter at different 
orders of chiral effective field theory and for different choices of the resolution scale. We discuss 
the convergence pattern of the predictions and their cutoff dependence. We explore to which 
extent fully polarized neutron matter behaves (nearly) like a free Fermi gas. We also consider 
the more general case of partial polarization in neutron matter as well as the presence of a small 
proton fraction. In other words, in our calculations, we vary both spin and isospin asymmetries. 
Confirming the findings of other microscopic calculations performed with different approaches, we 
report no evidence for a transition to a polarized phase of neutron matter. 
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I. INTRODUCTION 


The equation of state (EoS) of highly neutron-rich matter is a topic of current interest 
because of its many applications ranging from the physics of rare isotopes to the properties 
of neutron stars. In spite of recent and fast-growing effort, the density dependence of the 
symmetry energy, which plays a chief role for the understanding of those systems, is not 
sufficiently constrained and, at the same time, theoretical predictions show considerable 
model dependence. 

Polarization properties of neutron/nuclear matter have been studied extensively with a 
variety of theoretical methods [1-25], often with contradictory conclusions. In the study 
in Ref. [21], for instance, the possibility of phase transitions into spin ordered states of 
symmetric nuclear matter was explored based on the Gogny interaction [5] and the Fermi 
liquid formalism. In that paper, the appearance of an antiferromagnetic state (with opposite 
spins for neutrons and protons) was predicted, whereas the transition to a ferromagnetic 
state was not indicated. This is in contrast to predictions based on Skyrme forces [23] . 

The properties of polarized neutron matter (NM) have gathered much attention lately, 
in conjunction with the issue of ferromagnetic instabilities together with the possibility of 
strong magnetic fields in the interior of rotating neutron stars. The presence of polarization 
would impact neutrino cross sections and luminosities, resulting into a very different scenario 
for neutron star cooling. 

There are also other, equally important, motivations to undertake studies of polarized 
matter. In Ref. [26], for instance, we focussed on the spin degrees of freedom of symmetric 
nuclear matter (SNM), having in mind a terrestrial scenario as a possible “laboratory”. We 
payed particular attention to the spin-dependent symmetry potential, namely the gradient 
between the single-nucleon potentials for upward and downward polarized nucleons in SNM. 
The interest around this quantity arises because of its natural interpretation as a spin de¬ 
pendent nuclear optical potential, defined in perfect formal analogy to the Lane potential 
ra for the isospin degree of freedom in isospin-asymmetric nuclear matter (lANM). 

Whether one is interested in rapidly rotating pulsars or, more conventional, laboratory 
nuclear physics, it is important to consider both spin and isospin asymmetries. First, neu¬ 
tron star matter contains a non-negligible proton fraction. Concerning laboratory nuclear 
physics, one way to access information related to the spin dependence of the nuclear inter- 
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action in nuclear matter is the study of collective modes such as giant resonances. Because 
a spin unsaturated system is usually also isospin asymmetric, both degrees of freedom need 
to be taken into account. For those reasons, in previous calculations [2H], we extended our 
predictions |26l EH] to include matter with different concentrations of neutrons and protons 
where each nucleon species can have dehnite spin polarization. Our framework was based on 
the Dirac-Brueckner-Hartree-Fock (DBHF) approach to nuclear matter together with a real¬ 
istic meson-theoretic potential. Our hndings did not show evidence of a phase transition to a 
ferromagnetic (FM) or antiferromagnetic (AFM) state. This conclusion appears to be shared 
by predictions of all microscopic models, such as those based on conventional Brueckner- 
Hartree-Fock theory [T6|. On the other hand, calculations based on various parametrizations 
of Skyrme forces result in different conclusions. For instance, with the SLy4 and SLy5 forces 
and the Fermi liquid formalism a phase transition to the AFM state is predicted in asym¬ 
metric matter at a critical density equal to about 2-3 times normal density [23|. Qualitative 
disagreement is also encountered with other non-microscopic approaches such as relativistic 
Hartree-Fock models based on effective meson-nucleon Lagrangians. For instance, in Ref. |H| 
it was reported that the onset of a ferromagnetic transition in neutron matter, and its crit¬ 
ical density, are crucially determined by the inclusion of isovector mesons and the nature of 
their couplings. 

The brief review given above summarizes many useful and valid calculations. However, 
the problem common to all of them, including microscopic approaches, is that it is essentially 
impossible to estimate, in a statistically meaningful way, the uncertainties associated with 
a particular prediction, or to quantify the error related to the approximations applied in a 
particular model. 

Effective held theories (EFT) have shown the way out of this problem. Chiral effective 
held theory is a low-energy realization of QCD [30l|3T| which hts unresolved nuclear dynamics 
at short distances to the properties of two- and few-nucleon systems. Together with a power 
counting, chiral EFT provides a framework where two and few-nucleon forces are generated 
on an equal footing in a systematic and controlled hierarchy. 

Estimates of theoretical uncertainties [32| for calculations of the equation of state of 
nuclear and neutron matter have largely focused on varying the low-energy constants and 
resolution scale at which nuclear dynamics are probed [33VIHH] . In a recent work [301, we layed 
the foundations for order-by-order calculations of nuclear many-body systems by presenting 
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consistent NLO and N^LO chiral nnclear forces whose relevant short-range three-nncleon 
forces (3NF) are ht to A = 3 binding energies and the lifetime of the triton [32] • We 
then assessed the accuracy with which inhnite nuclear and neutron matter properties and 
the isospin asymmetry energy can be predicted from order-by-order calculations in chiral 
effective held theory. In this paper, we apply the same philosophy to study the equation of 
state of polarized neutron matter. 

Based on the literature mentioned above, a phase transition to a polarized phase (at 
least up to normal densities) seems unlikely, although the validity of such conclusion must 
be assessed in the context of EFT errors. Furthermore, polarized neutron matter is a 
very interesting system for several reasons. Because of the large neutron-neutron scattering 
length, NM displays behaviors similar to those of a unitary Fermi gas. In fact, up to nearly 
normal density, (unpolarized) neutron matter is found to display the behavior of an F-wave 
superhuid mnn. The possibility of simulating low-density NM with ultracold atoms near 
a Feshbach resonance |13] has also been discussed. When the system is totally polarized, 
it has been observed to behave like a weakly interacting Fermi gas |11]. Here, we wish to 
explore to which extent and up to which densities we are in agreement with such conclusions, 
and how this and other observations depend on the chiral order and the resolution scale. 

In comparison with the calculations of Ref. [H] (where 3NFs and 4NFs up to N^LO were 
included), our present work contains the following novelties: 

• We consider both cutoff dependence and truncation error for the purpose of uncer¬ 
tainty quantihcation of chiral EFT. Although incomplete in the 3NF at N^LO, our 
calculations are a substantial step in that direction. We note, further, that the contri¬ 
bution from the 3NF at N^LO was found to be very small in neutron matter for the 
potentials in our perview [321, about -0.5 MeV at normal density. Here, we consider 
neutron matter or highly neutron-rich matter. 

• For the hrst time, we present results for both spin and isospin asymmetries within the 
framework of chiral forces. As discussed in Section im these tools are necessary to 
assess, for instance, the sensitivity of the results (particularly, the potential onset of a 
phase transition) to the presence of a proton fraction. 

This paper is organized as follows: In the next section, we present the formal aspects 
of the self-consistent calculation of the energy per particle, in general, applicable to inhnite 
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matter with any degree of isospin and spin asymmetry. We also describe our approach to 
two- and three-body chiral forces. We provide expressions for the in-medium effective three- 
body force suitable for the most general case of different proton and neutron concentrations 
where each species can be polarized to a different degree. To the best of our knowledge, this 
has not been reported before in the literature within the framework of chiral forces. Results 
for polarized and partially polarized NM, as well as for polarized netron-rich matter in the 
presence of a small proton fraction, are discussed in Section III. Conclusions and future 
plans are summarized in Section IV. 


II. FORMALISM 


A. General aspects 


In a spin-polarized and isospin asymmetric system with hxed total density, p, the partial 
densities of each species are 

Pn Pnu T Pnd i Pp Ppu T Ppd i P Pn T Pp i (1) 

where u and d refer to up and down spin-polarizations, respectively, of protons (p) or neu¬ 
trons (n). The isospin and spin asymmetries, a, and fdp, are dehned in a natural way: 

Pn ~ Pp n Pnu ~ Pnd n Ppu ~ Ppd 
Q. , Pn , Pp . (zj 

P Pn Pp 

The density of each individual component can be related to the total density by 

Pnu (iT/^n) (iT*^) ^ t Pnd (1 /^n)(lT*^)^ i Ppu (^T'/3p)(l ^ ’ Ppd (1 /^p)(l 

(3) 

where each partial density is related to the corresponding Fermi momentum through pra 
=(A;]T)^/(67r^). The average Fermi momentum and the total density are related in the usual 
way as p = {2k%)/. 

The single-particle potential of a nucleon in a particular ra state, Ura, is the solution of 
a set of four coupled equations. 


Unu Unu,nu T Unu,nd T Unu,pu T Unu,pd 
Und Und,nu T Und,nd T Und,pu T Und,pd 


“1“ Ujiu^pd 

(4) 


(5) 
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Upu Upu^nu “1“ Upu^nd Upu^pu H“ Upu^pd 
^pd Upd^nu Upd^nd “ 1 “ Upd^pu “ 1 “ Upd^pd ; 


( 6 ) 

(7) 


to be solved self-consistently along with the effective interaction, the G-matrix. (The latter 
will be discussed in the next two subsections.) In the above equations, each Uraya' term 
on the right-hand side contains the appropriate (spin and isospin dependent) part of the 
interaction. Gray a'- More specifically. 


Uraya'(k) = ^ < TCT, T V |G(fc, (?) |r(T, t'ct' >, (8) 

where the summation indicates integration over the Fermi seas of protons and neutrons with 
spin-up and spin-down, and 


< r(j, rV'|G(fc, ^|(Tr, a'r' > 


Y. I < + a') > n < ^r; ^r'|T(r + r) > 1^ 

L,L',S,J,M,Ml,T ^ ^ z z 

X < LML,S{a + a')\JM >< L'Ml, S{a + a')\JM > 
Xt^'-^Yl,^^Xkrel)YLMpCKel) < LSJ\G{krel,Karn^^^^^ . (9) 


The G-matrix which appears in the formulas above is constructed from the two-nucleon 
potential and the effective density-dependent 3NF as explained later. 

The need to separate the interaction by spin components brings along angular depen¬ 
dence, with the result that the single-particle potential depends also on the direction of the 
momentum, although such dependence was found to be weak [28]. The G-matrix equation 
is solved using partial wave decomposition and the matrix elements are then summed as in 
Eq. (|^ to provide the new matrix elements in the representation needed for Eq. ([^, namely 
with spin and isospin components explicitely projected out. Furthermore, the scattering 
equation is solved using relative and center-of-mass coordinates, krei and Kc.m., since the 
former is a natural coordinate for the evaluation of the nuclear potential. Those are then 
easily related to the momenta of the two particles, k and q, in order to perform the inte¬ 
gration indicated in Eq. ([^. Notice that solving the G-matrix equation requires knowledge 
of the single-particle potential, which in turn requires knowledge of the effective interac¬ 
tion. Hence, Eqs.( |4]|^ together with the G-matrix equation constitute a rather lengthy 
self-consistency problem, the solution of which yields the single-nucleon potentials in each 
ra channel. 
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The kernel of the G-matrix equation contains the Pauli operator for scattering of two 
particles with two different Fermi momenta, and which is defined in analogy with 
the one for isospin-asymmetric matter |45j . 


Qr 


i{k,q,ky,k 


T <J 

F 


1 if p > ky and q > kp"^' 
0 otherwise. 


( 10 ) 


The Pauli operator is expressed in terms of krei and Kc.m. and angle-averaged in the usual 
way. We then proceed with the calculation of the energy per nucleon in the particle-particle 
ladder approximation, namely the leading-order contribution in the hole-line expansion. 
(See Ref. [ID] and references therein for a discussion of the uncertainty associated with this 
approximation.) 

Once a self-consistent solution for Eqs. (|4]j^ has been obtained, the average potential 
energy for a given rex component can be calculated. A final average over all ra components 
provides, along with the kinetic energy iPra, the average energy per particle in spin-polarized 
isospin-asymmetric nuclear matter. Specifically, 


f = 4 E E E + \ur,(k )). (11) 

^ ^ a=u,dT=n,pk<ky ^ 

where E/A is a function of p, a, /3„, and /3p, with q;= 1 in the present case. All calculations 
are conducted including values of the total angular momentum J from 0 to 15. 


B. Chiral two-body potentials 

In this section we discuss in some detail the features of the nucleon-nucleon (NN) po¬ 
tentials we use for these calculations. 

All low-momentum interactions are limited in calculations of the EoS to densities where 
the characteristic momentum scale (on the order of the Fermi momentum) is below the scale 
set by the momentum-space cutoff A in the NN potential regulating function, which for 
chiral NN forces typically has the form; 

f{p',p) = exp[-(p7A)2- - ip/Af-] , (12) 

where A < 500 MeV is associated with the onset of favorable perturbative properties [371 EH] ■ 
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NLO 

A (MeV) 

n 

Cl 

C3 

C4 


450 

2 





500 

2 





600 

2 




N^LO 

A (MeV) 

n 

Cl 

C3 

C4 


450 

3 

-0.81 

-3.40 

3.40 


500 

3 

-0.81 

-3.40 

3.40 


600 

3 

-0.81 

-3.40 

3.40 

N^LO 

A (MeV) 

n 

Cl 

C3 

C4 


450 

3 

-0.81 

-3.40 

3.40 


500 

2 

-0.81 

-3.20 

5.40 


600 

2 

-0.81 

-3.20 

5.40 


TABLE I: Values of n and low-energy constants of the dimension-two ttN 
each order and for each type of cutoff in the regulator function given in Eq. 
appears at NLO. The low-energy constants are given in units of GeV“^. 


Lagrangian, 01 ^ 3 ^ 4 , at 


(12). None of the Cj’s 


Although designed to reproduce similar NN scattering phase shifts, NN potentials with 
different regulator functions will yield different predictions in the nuclear many-body prob¬ 
lem due to their different off-shell behavior. On the other hand, appropriate re-adjustment 
of the low-energy constants that appear in the nuclear many-body forces is expected to 
reduce the dependence on the regulator function EH. 

In the present investigation we consider NN potentials at order (g/A^)^, (g/A^)^ and 
(g/A^)*^ in the chiral power counting, where g denotes the small scale set by external nucleon 
momenta or the pion mass and A^ is the chiral symmetry breaking scale. Chiral NN 
potentials at NLO and N^LO, corresponding to (g/A^)^ and (g/A^)^, have been constructed 
previously in Ref. 071 for cutoffs ranging from A = 450 MeV to about 800 MeV. With 
varying chiral order and cutoff scale, the low-energy constants in the two-nucleon sector are 
refitted to elastic NN scattering phase shihs and properties of the deuteron. The low-energy 
constants 0 ^ 3,4 associated with the ttttNN contact couplings of the chiral Lagrangian 
are given in Table We note that the c* can be extracted from irN ot NN scattering 
data. The potentials we use here [IHl SB] follow the second path. At N^LO, taking the 


















range determined in analyses of elastic 'kN scattering as a starting point, values were chosen 
to best reproduce NN data at that order. At N^LO, high-precision required a stronger 
adjustment of C 4 depending on the regulator function and cutoff. The htting procedure is 
discussed in Ref. |16], where it is noted that the larger value for C 4 has, overall, a very small 
impact but lowers the ^^2 phase shift for a better agreement with the phase shift analysis. 

In Ref. IH] , it was found that the two-body scattering phase shifts can be described well 
at NLO up to a laboratory energy of about 100 MeV, while the N^LO potential hts the data 
up to 200 MeV. Interestingly, in the latter case the y^/datum was found to be essentially 
cutoff independent for variations of A between 450 and approximately 800 MeV. Finally, we 
also use NN potentials constructed at next-to-next-to-next-to-leading order (N^LO) [461148] . 
with low-energy constants 04,34 displayed in Table 

Although N^LO calculations can achieve sufficient accuracy in selected partial wave chan¬ 
nels up to i?iab = 200 MeV, only the N^LO interactions achieve the level of high-precision 
potentials, characterized by a y^/datum ~ 1 . 

At the two-body level, each time the chiral order is increased, the NN contact terms 
and/or the two-pion-exchange contributions proportional to the low-energy constants 04 , 3,4 
are rehtted. We recall that at N^LO no new NN contact terms are generated, and therefore 
improved cutoff independence in the NN phase shifts [40] is due to changes in the two- 
pion-exchange contributions. At N^LO, subleading titiNN vertices enter into the chiral NN 
potential. These terms encode the important physics of correlated two-pion-exchange and 
the excitation of intermediate A(1232) isobar states. Therefore, at this order it is possible 
to obtain a realistic description of the NN interaction at intermediate range, traditionally 
generated through the exchange of a hctitious a meson of medium mass. At N^LO in the 
chiral power counting, 15 additional NN contact terms (bringing the total number to 24 at 
N^LO) result in a much improved description of NN scattering phase shifts. 

C. The three-nucleon force 

The leading three-nucleon force makes its appearance at third order in the chiral power 
counting and contains three contributions: the long-range two-pion-exchange part with 
ttttNN vertex proportional to the low-energy constants 04 , 03 , 04 , the medium-range one-pion 
exchange diagram proportional to the low-energy constant od, and hnally the short-range 
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FIG. 1: Diagrams for the chiral three-nucleon interaction at N^LO. In neutron matter, only diagram 
(a) contributes. 

contact term proportional to c^;. The corresponding diagrams are shown in Fig. labeled 
as (a), (b), (c), respectively. Diagrams (b) and (c) vanish in neutron matter, while all three 
terms contribute in symmetric nuclear matter [501151]. 

Although efforts are in progress to incorporate potentially important N^LO 3NF contri¬ 
butions [52H5l| both in the neutron and nuclear equations of state and the fitting of the 
relevant low-energy constants, the “N^LO” study reported in this paper is limited to the 
inclusion of the N^LO three-body force together with the N^LO two-body force, an approx¬ 
imation that is commonly used in the literature. The associated uncertainties for neutron 
matter have been investigated in Ref. m. 

To facilitate the inclusion of 3NFs in the particle-particle ladder calculation, we employ 
the density-dependent NN interaction derived in Refs. [l9l|50| from the N^LO chiral three- 
body force. This effective interaction is obtained by summing one particle line over the 
occupied states in the Fermi sea. Neglecting small contributions [SI] from terms depending 
on the center-of-mass momentum, the resulting NN interaction can be expressed in analyt¬ 
ical form with operator structures identical to those of free-space NN interactions, and are 
therefore included on the same footing as two-body forces. The small uncertainty associated 
with the use of these effective density-dependent 3NFs was discussed in Ref. [iO] . 

For the case of polarized isospin-asymmetric matter, the expressions from Ref. [SU] are to 
be extended to include four different Fermi momenta, namely those of upward(downward) 
polarized neutrons(protons), as described below. 

Using the notation established above to indicate the Fermi momenta of spin-up and spin- 
down neutrons or protons, the neutron and proton densities are given by pn = [{kp'^Y + 
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FIG. 2: Diagrams for the in-medium NN interactions corresponding to (i=l,...,6) given in 

the text. 


{kp^)^]/QTT'^ and Pp = [{k^)^ + (fc^'^)^]/67r^. 

Concerning kinematics, we consider elastic scattering process iVi(p) + A^ 2 (—p) —^ ^i(p + 
q) + N 2 {—p — q) in the center-of-mass frame. 

Following the notation of Ref. IHD], we can distingnish between six effective density- 
dependent NN interactions, denoted by diagram (1) to (6) in Fig. 2. They are: 

The Panli blocked pion-selfenergy (diagram (1)): 


■Yarned,! 9a ^ ^ (Tj • Q 0*2 ‘ 9 ^2 

^ NN r\ £ A "^2 / 2 I 2^2 

2/4 {mi + q^y 


( 2 cim| + C3(/^)(Pp + p„) , 


(13) 


The Panli blocked vertex correction (diagram (2)): 




med,2 

NN 


9a 


rD ■ T2 


ai-q(T2-q 
ml + 


16vr2/^ 

X g'(P+(p) + 2P+(p) + P+(p)) + 4P+(p) 


dcim^ [p+(p) + P+(p)] - (c3 -f C 4 ) 


+ 4c4 


2A{pp + pn) -mlVAp) 


9\ + 


327r2/4 
X 


ml + q^ 


- 4cim^ Pq (p) + Pi (p) + (C 4 - C 3 ) 


q (ro (P) + 2Pi (p) + P 3 (p)) + 4 P 2 (p)J + 4c4[27r {pn - Pp) + m^Pg (p) 

1 


9a 




:(r - t- 2 )*(<?! - ^ 2 ) ■{pxq) 


ml -|- 4p2 — g 2 
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(14) 


x|4cim^[ro(p)+ 1^ (p)J +C3(4p -q ) [Fq (p) + 21^ (p) + Tg (p) 

The last contribution, proportional to (rj^ — tI), leads to spin-singlet and spin-triplet mixing 
in the medium. It has been Fierz-transformed to bring it into the form of the anti-symmetric 
spin-orbit operator i{ai—a 2 )-{pxq). Terms which break rotational invariance in momentum- 
space due to the spin-polarization of the nuclear medium in 2 ;-direction have been discarded. 
Next, we give the expression for the Pauli blocked two-pion exchange (diagram (3)): 


4", 


medjS 

NN 


9a 


327r2/^ 


- 12cim^ 2r+(p) - {2ml + ?^)G'([(p,g) 


-3C3 


8vr^(p, + Pn) - 4(2m^ + q^)T^{p) - 2g^r+(p) + (2m^ + q^G^P, <D 


2-p+ / 


2'\2 yO+ / 


+4c 4 Ti ■ T 2 ((Ji ■ (J 2 - <?1 ■ g(T 2 ■ q)G 2 (p, g) 

-( 3 c 3 + Cifi ■ f2)i{ai + (T2) ■ {qxp) 2r+(p) 2r+(p) - {2ml + q^) 

X {Go {p, q) + 2G| (p, g))] - 12cim^ z(ai + (J 2 ) - {qxp) (p, g) + 2Gt (p, g) 
+4 c4 Ti ■ r 2 (Ti ■ (g X p) 0^2 ■ (g x p) (p, g) + 4G^ (p, g) + AG^ (p, g) 


9a 


+C 3 


(G + t- 2 )<^ 4cim^ 2 ro (p) - ( 2 m^ + q )Go (p, g) 


87r^(pp - pn) - 4{2ml + g2)ro (p) - 2q^V^ (p) + {2ml + q^fG^ (p, g) 

-4c 4 (ai ■ (T 2 g^ - ■ g <?2 ■ g )G 2 (p, g) 

-F(c 3 -1- C4) i{ai + (J 2 ) ■ {qxp) 2ro (p) -h 2rh (p) - {2ml + g^) 

X (Gq (p, g) + 2G^ (p, g))] + 4cim^ i{ai + ^ 2 ) - {qxp) [Gq (p, g) + 2G^ (p, g) 
-4 c 4 ai ■ (g X p) (J 2 ■ (g x p) [Gq (p, g) + 4Gh (p, g) + 4Gj (p, g) 


The loop functions r^(p) and Gf{p, q) with a superscript -|- or — are given by: 


rf (P) = \ kpu) + ri(P, kpd) 


i 2 ^j(.P^ ^nu) T ^^(p, knd) 


Gf{p, ?) = 2 [GjiP^ 9, kpu) + Gj{p, g, kpd) 


± 


2 L 


Gj{p, g, knu) + Gj{p, g, knd) 


(15) 

(16) 
(17) 


where Tj{p, kf) and G(p, g, kf) are dehned in Eqs. (13-16) and Eqs. (18-22) of Ref. 

Now we present the contributions from the Ivr-exchange 3NF proportional to the low- 
energy constant Cd- The vertex correction to Ivr-exchange linear in proton and neutron 
densities is (diagram (4)): 




med,4 


NN 


9aCd 

16^ 


2ri ■ T2{Pp + Pn) + (Ti r|)(Pp - Pn) 


ai- q (72- q 
ml + q^ 


(18) 
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Pauli-blocking, diagram (5), contributes: 




medjS 

NN 


QaCd 




Ti ■ r 2 


/ q 

2(Ti ■ (72 r+ (p) + f ai ■ (72 (2p^ - y■ g <?2 ■ g 


■ (<? X :p)^ 2 ■ (g X p)) r+(p) + 2r+(p) +r+(p) 


+ 127 r^(pp + Pn) - 

- (^1 + T-l) 12^1 ■ <?2 r2 (p) + 


QaCd / 3 3^ 


647r2/^A^ 


(7i ■ a2(2p^ - y) + <?1 ■ g(72 ■ g 


^ - 4 *^1 ■ ^ ^2 ■ (g X P) 

\ qz / qz 


+47r (p„ - pp) + 2m^ro (p) 


ro (p) + 2ri (p) + r3 (p) 


(19) 


The contribution from the contact 3NF proportional to the low-energy constant ce is 
(diagram (6)): 


T7med,6 _ 

^NN 


3ce 


- 2{pp + Pn) + (Pp - Pn)(ri^ + T^) 


( 20 ) 


Partial wave matrix elements with J > 1 of the antisymmetric spin-orbit term, which 
occur in Eq. (14), mix spin-singlet and spin-triplet states and these can be calculated for 
on-shell kinematics in the center-of-mass frame as: 


(J0J|*((7 i-(72 )-(p X q)F{q^)\JlJ) = (J1 J|*((7i- 0 ^ 2 ) ■ (p X q)F{p\ q^)\J0J) 


Phj + 1) p 

2J + 1 J-i 


dzp^F{p^, 2p^(l - z)) Pj-i{z) - Pj+i{z) 


( 21 ) 


However, because of the small size of this contribution, particularly for small proton frac¬ 
tions, we neglect this term in the present calculations. 


III. RESULTS AND DISCUSSION 

We show in Fig. [^the energy per particle in fully polarized neutron matter as a function 
of density. The yellow and red bands represent the predictions of complete calculations 
at second and third order, respectively, of chiral effective field theory, while the blue band 
shows the predictions obtained with the exploratory N^LO calculation as described above. 
For each band, the width is obtained by changing the cutoff between 450 MeV and 600 MeV. 

At N^LO and N^LO, cutoff dependence is generally moderate up to saturation density. 
At NLO, the cutoff dependence is practically negligible throughout. In unpolarized neutron 


13 



















matter, on the other hand, the largest cutoff dependence was seen at NLO [iQ]. This suggests 
that, in unpolarized NM, the larger cutoff sensitivity at NLO is mostly due to singlet states, 
particularly ^S'o, which are absent from the polarized system. At the same time, 3NFs do 
not appear at NLO, implying that most of the cutoff dependence in polarized NM at N^LO 
and N^LO is caused by the 3NF contributions. 

Clearly, the variations associated with changing the cutoff are not a good indicator of 
the uncertainty at a given order of chiral effective held theory, as the results from one order 
to the other do not overlap. Furthermore, the predictions do not show a good convergence 
pattern, although some indication of slow convergence can be seen when moving from N^LO 
to our N^LO calculation. 

As can be concluded from Table the predictions from the N^LO calculation are close 
to the free Fermi gas energy, at least up to saturation densities. Our results with the N^LO 
(A=500 MeV) potential are in good agreement with those from Ref. [H] using the same 
potential as well as three- and four-nucleon forces at N^LO. With regard to the similarity 
with the free Fermi gas, it is interesting to include some additional considerations. As 
mentioned in the Introduction, many-fermion systems with large scattering lengths offer the 
opportunity to model low-density neutron matter. In the unitary limit (that is, when the 
system can support a bound state at zero energy), the scattering length approaches inhnity. 
The system then becomes scale-independent and the ground-state energy is determined by a 
single universal parameter, known as the Bertsch parameter, The latter is dehned as the 
ratio of the energy per particle of the unitary gas to that of the free Fermi gas. In Ref. [53] , 
using a simple ansatz for the interaction, it is shown that ^ increases from approximately 
0.5 to 1.0 as the spin asymmetry of neutron matter, is increased from 0 (unpolarized) 
to 1 (fully polarized). 

In Fig. 1^ for our N^LO calculation, we compare predictions (along with their cutoff 
variations) of the energy per neutron in: unpolarized NM (green band), partially polarized 
NM (pink band), and fully polarized NM (blue band). For the partially polarized case, the 
value of /3n (see Eq. (2)) is equal to 0.5, corresponding to 75% of the neutrons being polarized 
in one direction and 25% in the opposite direction, see Eqs. (3-4). Clearly, a lesser degree of 
spin asymmetry (as compared to the ferromagnetic case) yields considerably less repulsion. 
There is dehnitely no sign of a phase transition, particularly to a ferromagnetic state, nor 
an indication that such transition may occurr at higher densities. This is consistent with 
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Density (fm 

A (MeV) 

Effg/E 

0.15 

450 

0.95 


500 

0.92 


600 

0.95 

0.17 

450 

0.95 


500 

0.91 


600 

0.93 


TABLE II: Ratio of the energy per particle of a free Fermi gas to the energy per particle of 
polarized neutron matter around saturation density at N^LO (as described in the text) and for 
different values of the cutoff. 

what we observed earlier [29] with meson-theoretic interactions. 

As a baseline comparison, we also include, for the unpolarized case, predictions based on 
a different approach, shown by the black dotted line in Fig. These are taken from Ref. [SH] 
and are based on the Argonne uis two-nucleon interaction plus the Urbana IX three body- 
force, using variational methods. The predictions are overall in reasonable agreement with 
our green band, although those from Ref. [56] show more repulsion as compared to the softer 
chiral interactions. 

Most typically, models which do predict spin instability of neutron matter hnd the phase 
transition to occurr at densities a few times normal density. Such high densities are out¬ 
side the domain of chiral perturbation theory. With some effective forces, though, it was 
found UTl that a small fraction of protons can signihcantly reduce the onset of the threshold 
density for a phase transition to a spin-polarized state of neutron-rich matter. We explored 
this scenario by adding a small fraction of protons to fully polarized or unpolarized neutrons. 
From Eqs. [I]j^ a proton fraction of 10% is obtained with a;=0.8. The results are displayed 
in Fig. [^ where a crossing of the bands labeled with “0.8, 1.0” and “0.8, 0.0”, respectively, 
would indicate a phase transition. Thus we conclude that such transition is not predicted 
with chiral forces. By extrapolation, a transition to a polarized state would also appear very 
unlikely at higher densities. 
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FIG. 3: (color online) Energy per neutron in fully polarized neutron matter as a function of 
density. The yellow and red bands represent the uncertainities due to cutoff variations obtained in 
the complete calculations at NLO and N^LO, respectively. The blue band is the result of the same 
cutoff variations applied to our exploratory N^LO calculation, see text for details. The dotted 
curve shows the energy of the free Fermi gas. 

IV. CONCLUSIONS AND OUTLOOK 

We have calculated the equation of state of (fully and partially) polarized neutron-rich 
matter. We performed complete calculations at second and third order of chiral effective 
held theory and calculations employing the N^LO 2NF plus the leading 3NF. Results with 
both spin and isospin asymmetries are presented for the hrst time with chiral forces. 

In all calculations, the cutoff dependence is moderate and dehnitely underestimates the 
uncertainty of each order. Concerning the latter, we do not see a satisfactory convergence 
pattern. The missing 3NFs are most likely not the main cause of uncertainty at N^LO, since 
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p (fm-3) 


FIG. 4: (color online) Energy per neutron in pure neutron matter as a function of density at 
N^LO. From lowest to highest curve: unpolarized NM; partially polarized NM, with /3„=0.5; fully 
polarized NM (/3„=1). The width of each band shows the uncertainty from varying the cutoff 
between 450 and 600 MeV. The black dotted line shows the predictions for the equation of state 
of unpolarized neutron matter from Ref. |56j . 

Ref. [H] has demonstrated that large cancelations take place between the 27r-exchange 3NF 
and the vr-ring 3NF at N^LO, while other 3NF contributions are very small (about 0.1-0.2 
MeV). Clearly a calculation at N'^LO is absolutely necessary to get a realistic indication of 
the EFT error at N^LO. Such effort is in progress. If such calculation displays a reason¬ 
able convergence pattern, it will be strong evidence that polarized neutron matter, indeed, 
behaves nearly like a free Fermi gas, at least up to normal densities. 

In our N^LO calculation, the energies of the unpolarized system at normal density are 
close to 16 MeV for all cutoffs, whereas those in the polarized case are approximately 60 MeV. 
Thus, even in the presence of the large uncertainties discussed above, a phase transition to 
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100 



FIG. 5: (color online) Energy per nucleon in neutron-rich matter as a function of density at N^LO 
and different conditions of isospin and spin polarization. The (brighter blue) band labeled as “0.8, 
1.0” displays the results for neutron-rich matter with a proton fraction equal to 10% (q:=0.8) 
and fully polarized neutron (/1„=1.0). The (brighter green) band labeled as “0.8, 0.0” refers to 
neutron-rich matter with the same proton fraction and no polarization (/3„=0.0). The protons are 
unpolarized. For comparison, we also include the bands (darker blue and darker green) already 
shown in the previous figure, which refer to pure neutron matter (a=l) with fully polarized (/3n=l) 
or unpolarized (/3„=0) neutrons. The bands are obtained varying the cutoff between 450 and 600 
MeV. 

a ferromagnetic state can be excluded. This conclusion remaind valid in the presence of a 
small proton fraction. 
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